NASA TM X- 53288 


June 28, 1965 


oo 

oo 


I 

X 

s 

H 

CD 

< 


DETERMINATION OF AERODYNAMIC FORCE 

AND HEAT TRANSFER PROPERTIES FOR A CONCAVE 

HEMISPHERICAL SURFACE IN FREE MOLECULAR FLOW 


N65-30184 


(ACCESSION NUMBER) 

4ter 

(THRU) 

/ 

(PAGES) 

(CODE) 

M. 

(NASA CR OR TMX OR AD NUMBER) 

(CATEGORY) 


NASA TECHNICAL 
MEMORANDUM 


by CLARENCE RAY WIMBERLY 
Aero-Astrodynamics Laboratory 


NASA 

George C. Marshall 
Space Flight Center , 

Huntsville , Alabama 


GPO PRICE $ 
CFSTI PRICE(S) $ 


Hard copy (HC) 
Microfiche (MF) 

ff 653 July 65 


& cL 


TECHNICAL MEMORANDUM X-53288 


DETERMINATION OF AERODYNAMIC FORCE AND HEAT TRANSFER PROPERTIES FOR A 
CONCAVE HEMISPHERICAL SURFACE IN FREE MOLECULE FLOW 


By 

Clarence Ray Wimberly 


ABSTRACT 


The methods and equations are presented for the aerodynamic force 
and heat transfer properties of a concave hemispherical surface at an 
arbitrary angle of attack in free molecule flow. The effect of multi- 
ple reflections is taken into account by assuming a cosine distribution 
after each collision. It is assumed that the emission of molecules from 
the surface is purely diffuse, and that the reflecting molecules are 
perfectly accommodated to the surface conditions. The resulting equations 
are not in closed form, but have been solved numerically for some typical 
values of the molecular speed ratio and surface-to-ambient- tempera ture 
ratio. It is found that under the above conditions the heat transfer 
characteristics are the same as those of a convex hemisphere. The drag 
is increased slightly over its counterpart, while the total lift is 
shown to be very small and would be zero except for interref lected 
molecules . 
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DEFINITION OF SYMBOLS 


Symbol 

A 

A r 

a 

s 

C T 


D 

d 

Ei 

E r 

E w 

f 

F a 

k 

m 

Ni* 


N, 


w 


Definition 

surface area of hemisphere (t:D 2 /2) 
reference area (:n:D 2 /4) 
thermal accommodation coefficient 
aerodynamic drag coefficient 

aerodynamic lift coefficient 

specific heat at constant pressure 
heat transfer coefficient 
specific heat at constant volume 

relative surface velocity in the x-direction 
relative surface velocity in the y-direction 

relative surface velocity in the z-direction 
aerodynamic drag and surface diameter 
distance between two differential surface elements 
energy due to incident molecular flux 
energy due to reflected molecular flux 

energy due to molecules reflecting at the wall temperature 

Maxwells velocity distribution function 

angle factor 

Boltzmann* s cons tant 

molecular weight 

total molecular incident flux 

total reflected molecular flux at the wall temperature 


iv 


Symbol 

n 


P 




DEFINITION OF SYMBOLS (Continued) 

Definition 

molecular density 

total pressure on a differential surface element 

pressure on a surface element due to the frees tream 
incident and reflected molecular flux 

pressure in region I 
pressure in region II 

pressure due to incident frees tream molecular flux 

pressure due to reflected frees tream molecular flux 

pressure due to reflected frees tream molecular flux at 
the wall temperature 

local heat transfer to the surface per unit time and area 

average heat transfer to the surface 

gas constant 

radius of hemisphere 

ratio of the free stream velocity and the most probable 
molecular speed 

molecular speed ratio at wall temperature 
molecular temperature in frees tream 

wall temperature 


v 



Symbol 
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DEFINITION OF SYMBOLS (Continued) 


v 

a 

P 

7 

0 

v 

JT 


P 

a 


T. 

1 

T r 

02 


Definition 

surface frees tream velocity 

molecular velocity 

surface angle of attack 

angle of integration 

ratio of specific heats (C /C ) 
r v p v 

element angle of attack 
angle of integration 
constant (3.14159) 
mass density 

reflection coefficient due to shear stresses 

reflection coefficient due to pressures 

shear stress due to incident mass flux 

shear stress due to reflected mass flux 

direction angle for diffuse emission from element 1 

direction angle for diffuse emission from element 2 
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DETERMINATION OF AERODYNAMIC FORCE AND HEAT TRANSFER PROPERTIES FOR A 
CONCAVE HEMISPHERICAL SURFACE IN FREE MOLECULE FLOW 


SUMMARY 


3 o/S4 


An analysis of a concave hemispherical surface in free molecule flow 
is presented where the surface is oriented at an arbitrary angle of attack. 
General equations are formulated for the lift, drag, and heat transfer 
coefficients. The effect of multiple reflections is taken into account 
by assuming a completely diffuse reflection process where molecules are 
distributed equally in all directions after each reflection. The mole- 
cules are assumed to obey Maxwell's velocity distribution law before and 
after collisions, reflecting directionally in a cosine distribution. It 
is further assumed that the reflecting molecules are perfectly accommoda- 
ted to the surface conditions. The resulting equations are not in closed 
form, but have been solved numerically for some typical values of the 
molecular speed ratio and surface-to-ambient- temperature ratio. It is 
found that under the above conditions the heat transfer characteristics 
are the same as those of a convex hemisphere. The drag is increased 
slightly over its counterpart, while the total lift is shown to be very 
small and would be zero except for interref lected molecules. Maximum 
drag and heat transfer occur at zero angle of attack, decreasing to zero 
at ninety degrees angle of attack. The lift is a function only of inter- 
reflected molecules and is zero at zero and ninety degrees angle of attack, 
the maximum occurring near an angle of attack of forty-five degrees. 



I. INTRODUCTION 


With the advent of orbital satellites and space vehicles, analyses 
of the aerodynamic force and heat transfer properties of these bodies and 
vehicles, in the high altitude, low density flow regime have become 
necessary. The complexity of the body shapes, and the peculiarities of 
the mqlecular interactions with surfaces in the flow have introduced 
interesting aspects to analyses in this flow regime by imposing the 
necessity of considering various contributions of the momentum and energy 
imparted to the body by the impinging and re-emitting molecules. Actually, 
much work has been done in determining the aerodynamic properties of flat 
and convex bodies in free molecule flow with experimental verification, 
but little has been done in determining these properties for concave bodies 
that create a molecular multi-reflection process. This paper is con- 
cerned with the determination of the aerodynamic properties of a concave 
hemispherical surface at various angles of attack in this high altitude, 
free molecule flow regime. 

By definition, the free molecule flow regime is a gas region of 
extreme rarefaction where collisions between molecules are much less 
probable than are collisions between the molecules and a surface in the 
flow. Because of this molecular freedom, the gas does not behave as a 
continuous fluid, but rather performs in accordance with its complex 
molecular structure defined by the kinetic theory of ideal gases. 

Briefly, this theory is a statistical description of molecular motion, 
including the impact effect with other bodies in the flow. From this 
theory and Newton* s laws of motion, values for the force and heat transfer 
properties imparted to a concave hemispherical surface submerged in the 
flow have been determined and are presented here. The equations were 
derived by treating the flows of incident and re-emitted molecules 
separately, where the incident flow consists of the frees tream and reflected 
molecules. It is the inclusion of this additional incident flux created 
by the concavity of the surface that causes the interesting aspects of 
this inves t iga t ion , 
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II. PREVIOUS WORK 


The fundamental concept of the free molecule gas flow analysis per- 
formed here is based on textbooks by Jeans [1], Loeb [2], and Patterson 
[3], the paperback by Schaaf and Chambre [4], and the paper by Oppenheim 
[5]. The methods employed in the analysis of free molecule flow over 
non-convex surfaces are presented in the book by Jakob [6], and the 
papers by Chahine [7], Sparrow [8], and Pratt [9]. 

The fundamentals of particle dynamics and the kinetic theory of 
gases, including velocity distributions, pressures, heat conduction 
and interaction effects are fully discussed with derivations in the 
textbooks by Jeans [1] and Loeb [2]. Using these fundamentals, Patterson 
[3] and Schaaf and Chambre [4] derived the means for representing the 
force and heat transfer properties for a differential surface element 
and for a few flat and convex bodies. Patterson [3] limited his 
investigation to the drag and heat transfer properties of a cylinder 
in free molecule helium and discussed some effects of Knudsen number on 
heat transfer. Schaaf and Chambre [4] summarized the pertinent equa- 
tions in deriving the aerodynamic force and heat transfer properties 
for a flat plate, cylinder, and sphere. Also included were data for the 

drag coefficient, lif t-to-drag ratio, recovery factor and Stanton number 

/ 

for the flat plate at angles of attack, and for the drag coefficient, 
recovery factor and Stanton number for the cylinder, and sphere, each 
plotted as a function of the molecular speed ratio. Oppenheim [5] dis- 
cusses in detail the derivation and application of the convective heat 
transfer equations to bodies in free molecule flow. He included specific 
equations and data for the heat transfer properties to flat plates, 
cylinders , semicylinders , spheres, and hemispheres. 

Jakob [6], in Chapter 31, discusses diffuse radiation between two 
surfaces of arbitrary orientation. The results of his discussions and 
derivations produce an equation representing the fraction of energy 
emitted from one surface to another. Chahine [7], Sparrow [8], and 
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Pratt [9] used this in determining the force and heat transfer properties 
of concave surfaces. Chahine [7] determined the force and heat transfer 
properties for infinite concave cylinders and for concave spherical 
segments; however, he limited his investigations to the special case 
wherein the free -stream flow impinges directly on all parts of the con- 
cave surfaces. This is an especially strong restriction for hemispherical 
surfaces since, for this type of surface, the angle of attack must be 
zero. Sparrow [8] in his paper determined the local and average heat 
transfer rates, the adiabatic wall temperature, and the forces exerted 
on a concave cylindrical surface in free molecule flow at arbitrary 
angles of attack. Pratt [9] performed the same analysis as Chahine [7], 
obtaining slightly different results. He also restricted his investiga- 
tion for hemispherical surface to zero angle of attack. The investiga- 
tions by Chahine [7], Sparrow [8], and Pratt [9] all assumed that the 
surface was traveling at hypervelocities, thereby neglecting the molecular 
motion, and assumed that all molecular reflections were diffuse. 

The analysis presented here is for the force and heat transfer proper- 
ties of concave' hemispherical surface moving at hypervelocities in free 
molecule flow including the effects of angle of attack from zero to ninety 
degrees. Perfect thermal accommodation of molecules at the surface is 
assumed where the molecules are diffusely reflected. 


Ill . FUNDAMENTAL CONCEPTS 
A. Free Molecule Flow 

The free molecule flow regime may be defined mathematically 
as the gas region where the Knudsen number K n = A/L ^ 10 [3]. This 
corresponds to a minimum altitude in the earth's atmosphere of about 
75 miles for a characteristic length of one foot, based on the U. S. 
Standard Atmosphere, 1962 [10]. This is an arbitrary definition and 
expressed by some as not complete for defining this flow regime [4] . 
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The Knuds en number, however, has been accepted by most as the 
best criterion to date for establishing the boundaries of free molecule 
flow and is defined as the ratio of the mean free path of the molecules 
(A) and some characteristic length (L) . The mean free path is the 
average distance a molecule travels before colliding with another 
molecule, 

B. Maxwells Velocity Distribution 

Because of the character of free molecule flow the methods 
normally used in describing fluid motion in a continuum cannot be 
applied. Thus, in this region of extreme rarefaction, the molecular 
motion is necessarily described statistically because of the variation 
of molecular velocities in the flow. From the kinetic theory of 
gases [1,2], the number of molecules in a region of free molecule 
flow is a function not only of the volume of the region, but also of 
the velocity of each molecule in the region. This differential number 
of molecules is thus expressed as 

dN = f(x,v,t) dx 3 dv 3 , (1) 

where the function f(x,v,t) has the same characteristics as molecular 
density except in a six- dimensional phase space. By assuming that the 
gas region of investigation is homogeneous and steady, the function 
becomes independent of position and time. The derivation of the result- 
ing function is presented in many books on gas kinetics such as those by 
Jeans [1] and Loeb [2]. This function is 


f(v) 




( 2 ) 


and is known as Maxwell's velocity distribution function [1,2, 3, 4]. 

As will be shown, this function is highly useful in determining the 
force and heat transfer properties of surfaces in free molecule flow. 
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C. Reflection and Accommodation Coefficients 


Specification of molecule-surface interactions is necessary to 
determine the momentum and energy transfer to the surface. Fortunately, 
this can be done by simply knowing certain average parameters which 
characterize the interaction phenomena. These parameters are the well 
known reflection and accommodation coefficients [3,4]. Two parameters 
are used in defining the reflection process to take into account both 
the tangential and normal force components. These reflection coefficients 
are . 


o 



( 3 ) 


a 


f 



( 4 ) 


For a purely specular reflection process, a = a 1 =0, while for a com- 
pletely diffuse reflection process, o = a 1 = 1 . The thermal accommoda- 
tion coefficient is defined by 


dE. - dE 
i r 

dE. - dE 
i w 


( 5 ) 


and is used in defining energy exchanges with the surface. From this 
definition, it is easily seen that for perfect thermal accommodation at 
the surface, a = 1, while for the case of no energy exchange, a = 0. 
Experimentally determined values of 0 , 0 f and a for various typical 
materials and air-surface combinations are available [3,4], 
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D. Pressure and Shear Stress Calculations 


The determination of the pressure and shear stress on a dif- 
ferential surface element (Figure 1) caused by the incident freestream 
and reflected molecules, where the macroscopic velocity is in the xy-plane, 
can now be carried out by considering the time rate of change of momen- 
tum normal and parallel to the surface element. By use of equations (1) 
and (2), the pressure and shear stress caused by the incident freestream 
are obtained as follows: 


P. 

i 



m C 2 f(v) dC 
x x 


-oo -oo O 



dC 

z 


— — — [7 + f) 2 ] [1 + erf (q) ] + qe -T1 "j- (6) 

2j~n S 2 l 2 J 


and 


00 00 00 



-00 -00 o 


f(v) 


dC dC dC 
x y z 


pU^ cos e r,-f + erf(n)1 l ( 7 ) 

2 ■JT S l -> 


where 

T) = S sin 0 
S = U A/ 2RT 
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erf On) = 



dx 


C = U + v 

XXX 

C = U - v 

y y y 

c = v . 

z z 

The limits of integration include all possible magnitudes of the 
molecular velocities. The lower limit of dC^ is zero since no negative 
velocities in the x-direction will strike the front side of the surface 
element. The function f(v) is given by equation (2). The angle 0 is 
the element angle of attack. The net pressure and shear stress can be 
determined by using the definitions of the reflection coefficients, 
equations (3) and (4): 


P r = P. + P = (2 - a 1 ) P. 4- a 1 P 
fir l w 


T f 


T. 

1 


T 


r 


OT i , 


( 8 ) 

(9) 


where P w is the pressure exerted on the wall due to molecules reflect- 
ing at the wall temperature with no macroscopic velocity contribution, 
i.e., S = 0 and T = T w . With the help of equations (1) and (6), the 
equation for p w becomes [4] 

P - i m ■7" 2 jcRT N., (10) 

w 2 w x 
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where 


00 CO oo 



C f(v) dc dC dC dA 
x x y z 



4kF/ 2jt Te” +n/7 t)[1 + erf(T])]l 


l 


dA. 


( 11 ) 


By combining equations (6), (7), (8), (9), (10) and (11), the net pres- 
sure and shear stress from the frees tream incident flux of molecules on 
the front side of the surface element dA (Figure 1), including the contri- 
bution of reflecting molecules, become 



T f = 


np U 2 cos 0 


2 \TjT ! 


j'e’ 7 ' + n/jT T) [1 + erf(q)]j- . 


(13) 


E. Heat Transfer Calculations 

The fundamental concepts of heat transfer to a differential 
surface element in free molecule flow is carried out in detail by 
Patterson [3], Schaaf and Chambre [4], and Oppenheim [5], where the 
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contributions of energy imparted to the surface are the translational 
kinetic energy and internal molecular energy. The ‘total kinetic energy 
caused by translational motion of molecules that strike the front side 
of dA per unit time and per unit area is 


dE. . 



4 mC* f (v) C dC dC dC 
2 x x y z 


= pRT n/rt/ 2at |(S 2 + 2)e' J]S + -JV (S 2 + |) T)[l + erf (rtf ]j-. (14) 


From the kinetic theory of gases [1,2] and .the principle of equipartition 
of energy, the internal energy per molecule is jkT/2, where j is the 
internal degrees of freedom. The flux of internal energy per unit time 
and per unit area is given by 


dE. . = jyN, 
1,1 J 2 i 


(15) 


where 

i 



*L 

1 


and is defined by equation (11). From the definition of the 
accommodation coefficient, equation (5), the heat transfer to the 
surface element becomes 


q f = dEL - dE r = a(dE t - dEj , 


(16) 


I 
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where 


dE. = dE. + dE. , 
i i,t 1,1 


( 17 ) 


and dE^ is the energy transfer from the surface caused by molecules 
reflecting at the wall temperature with no macroscopic velocity contri- 
bution, i.e., S = 0 and T - T w . Inspection of equations (11) and (14), 
while considering the contribution of internal energy, reveals the 
equation 


dE 

w 


y Hh 1 

2(7 - 1) 


mRT 

w 


N . 


w 


(18) 


Inserting equations (14), (15), and (18) into equation (16) gives the 
heat transfer to the surface element (Figure 1) caused by freestream 
incident and reflected molecules on the front side of the surface per 
unit time and per unit area 


q f = apRT 




y+1 W 
'2( 7 - 1) T 


^ + n/jT t )[1 + erf(r))] 


. i a ”! 2 


(19) 


IV. EFFECTS OF CONCAVITY 

To determine the effect of the concavity of the surface, it is 
necessary to represent the fraction of reflecting molecules that re-impinge 
on the surface and the effect on the mass and energy transfer to the sur- 
face. This may be done by assuming a diffuse reflection process which 
is based on Lambert's cosine law of diffusion and is defined simply as 
a cosine distribution of radiation or mass flux to all directions of 
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space [6]. More precisely, from Figure 2, it can be seen that from the 
surface element dA 1? mass is diffusely emitted in the direction QM pro- 
portional to the cosine of the angle 0 between OM and ON, the normal to 
dA x * For concave surfaces, a fraction of molecules reflect or reemit, 
either diffusely or partially diffusely, from portions of the surface and 
impinge on other parts of the surface, thereby creating additional con- 
tributions to the forces and heat transfer properties of the body. This* 
molecular fraction of reimpinging molecules will be referred to through- 
out this paper as the angle factor. The angle factors for a concave 
hemispherical surface at an arbitrary angle of attack can be determined 
by considering two differential surface elements as shown in Figure 3. 
Since the diffuse interchange of therma.l radiation also applies to the 
diffuse interchange of mass [6,8], the general relationship for the 
angle factor derived by Jakob [6] for mass directly incident (normal) 
to the surface element can be used 


F = 
a 


cos 0T cos 0p dA 
ltd 2 


( 20 ) 


For concave hemispherical surfaces, it can be easily seen from Figures 4 
and 5 that 

„ nf d 

cos 0 ± = cos 0 2 = — 


and 


dA 2 = r 2 sin {3 dp dv 2 * 


Thus, for one reflection, 


F 

a 


it 




( 21 ) 
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As may be shown from Figure 4, v Q - Jt - 2a . The angle factor for one 
reflection can thus be expressed as 



where the surface angle of attack, a, is measured in radians. This quantity 
represents the fraction of energy caused by incident molecules from one reflec- 
tion that re-impinge on a differential surface element at any point on 
the hemispherical surface. For a greater number of collisions the angle 
factor becomes a fraction of the molecules that was previously reflected. 

The general expression for the angle factors for multi-reflections can thus 
be written 



(23) 


where j = 1,2,3, ... , representing the number of interreflections (one 
less than the number of collisions). 

These angle factors defined by equation (23) represent only the 
component of energy caused by molecules striking normal to the surface 
element. The tangential component is not given since, by definition of 
diffuse emission, reflecting molecules are distributed equally and in 
opposite directions, thereby cancelling any effect they may have. 

V. ANALYSIS OF A CONCAVE HEMISPHERICAL SURFACE 

Based on the fundamental concepts as presented in Section III and 
the concavity effects of Section IV, the force and heat transfer can 
be determined for a concave hemispherical surface at an arbitrary angle 
of attack. A completely diffuse reflection process is assumed with per- 
fect thermal accommodation; thus, a = a 1 =1, which is considered justi- 
fied by experimental results [4], 


13 


A, Force Coefficients 


The force coefficients can be determined by integrating the 
pressure contributions of a surface element over the hemispherical sur- 
face. This pressure includes that caused by the frees tream molecular 
flux as defined by equation (12), plus contributions of reflecting 
molecules. For clarity, the pressure distribution over the surface can 
be divided into two regimes (see Figure 4). The area of direct incidence 
is defined by 0 I v ^ v Q and is termed as region I, while the area 
receiving only reflected molecules is defined by v Q ^ v ^ rt and is 
termed as region II. The dividing line is a semi-circle and can be 
easily shown to be at the angle 

v G = it - 2a, (24) 

where the surface angle of attack, a, appears only in the xy-plane as 
shown in Figure 1. No generality is lost in this assumption since the 
surface is axisymmetric. The pressure at any point in region I can be 
written as 


P 


I 




(l/2) j "\ 


(25) 


where Pr = Pw> since a 1 = 1. The value of P w is defined by equation (10). 
It can be shown that the infinite series in equation (25) converges, and 
that the sum of the series converges to 

00 

^ (l/2) j_1 = 2. 

j=l 
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Equation (25) can thus be written as 


P T = P- + 2 4 - % ^2 rt RT N , 

lx ^ ft W 1 


( 26 ) 


where and P^ are defined by equations (11) and (12), respectively. 
Similarly, the pressure at any point in region II can be determined to 
be the same as equation (26) minus the frees tream contribution. The 
pressure on the surface in region II can thus be written 

P II = 2{ \ " N.. (27) 

The differential force on the surface in the direction of flow measured 
by the angle of attack a, can now be determined. The derivation is 
facilitated by reference to Figure 4. This force includes the contri- 
butions of both pressure and shear stresses plus the shading effect due 
to angle of attack. 


dD = (Pj sin 0 dAj - T f cos 0 dAj + Pjj sin 0 dA-^j) sin (3 (28) 

where Pj and Tf are defined by equations (26) and (13), respectively, 
and where P^ is defined by equation (27). Also required are the dif- 
ferential surface areas , which are shown to be 

dA^ = dA^j = r 2 d(3 sin (3 dv. (29) 

The element angle of attack is defined as 0 = n - CC - v. The total drag 
or force, in the direction of flow can be determined by integrating equa 
tion (28) over the hemispherical surface. 
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V 

at o 


■II 


o o 


V 2 <2 


% m 'J 2itRT N. 
it w i 


sin 0 r 2 sin 2 


(3 dvdp 



+ 


it it 



e r 2 sin 2 P dvdp 


o 


V 


o 


(30) 


As previously defined the values of Pf, N^, Tf and are shown to be 
a function of the element angle of attack, 9. The appropriate substitu- 
tion into equation (30) can thus be made for v, changing the limits of 
integration. Furthermore, by inserting equations (11), (12), (13), and 
(27) into the above equation, where o’ = a = 1, the drag equation becomes 


JT-CC 


D = P RT ^ r J [g(e) + h(e)] sin 0 de + 


pRTnr £ 


it-a 


f(0) cos 2 0 d0 


a 


a 


pRTnr 2 (l 


n-a 


(31) 


- a 


where 


g(e) = 


U7 


S sin 0 + 


1-/VT 


-S 2 s in 2 0 


(32) 
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h(9) = [1 + erf (S sin 0) ] 


\ + 


. 2 „ , S sin 0 I ~ r~ 

sin 0 + - — -\/ irT /T 

2 w 


(33) 


f(0) = e +v7 S sin 0 [1 + erf(S sin 0) ] . 


(34) 


Similarly, the lift, or force normal to the direction of flow can be 
determined. Utilizing Figure 4, the differential lift can be determined 
as 


dL = (Pj cos 0 dA^. + sin 0 dA^. + P cos 0 dA^) sin(3 (35) 


where Pj, Tj, and Pj-j retain the same definitions. A similar analogue 
can be applied for the lift as was applied for the drag, yielding the 
following equation 


v 

Jt O 

L = J T f sin 0 r 2 sin 2 p dvdp 
o o 


v 

Jt o 


If [v0-$ 


- ) m \l 2jtRT N. 

Jt/ w 1 


cos 0 r 2 sin 2 p dvdp 


o o 


Jt Jt 

cos 0 r 2 sin 2 p dvdp. 

o v„ 


(36) 
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By making the same substitutions as was previously done for the drag equa- 
tion the resulting lift equation becomes 


n-a 


jt-a 


L = 


pRTjtr" 


[g (©) + h(e)] cos e de - 


pRTjtr^ 


f(0~) sin 0 cos 0 d9 


a 


a 


rc-a 


pRTnr 2 


\ - ^ n/t w /t J f(e) cos e de. 


(37) 


•a 


The force coefficients C D and C L for the hemisphere can thus be written 


CL = 


D 




it -a 

i r 

YP J tg( 0 ) + Me)] sin 0 d 0 

a 


jt-a 


n-a 




f(e) cos ^ e de + -^s - ~j 


7TP (i - JtJt r sin e de 


a 


^ ' * j 

-a 


(38) 


n-a 


c t = T 

L i pU 2 A 

2 h r 


2S‘ j 

a 


r 


[g(e) + Me)] cos 0 d0 


jt-a 


jt-a 


>777^ 

a 


r 


f(0) sin 0 cos 0 d0- — 


a 


, n/t/T 
n/ w 


-a 


COS 0 d0 


(39) 
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B. Heat Transfer Coefficients 


The heat transfer properties considered here for the concave 
hemispherical surface is based on the general theory of frees tream con- 
vective heat transfer in free molecule flow presented in the paper by 
Oppenheim [6], plus the contributions due to reflected molecules. The 
heat transfer to a differential surface element due to the incident flux 
is defined by equation (19). The heat transfer due to reflected molecules 
is determined by using the angle factors of Section' IV where the energy 
impinging on a differential element is equal to molecular energy emitted 
from the rest of the surface times the appropriate angle factor. For 
example, if dE r leaves element dA ly then the incident amount that is 
transferred to dA 2 is F a dE r . The initial heat transfer caused by the 
incident freestream and reflected molecules is defined by equation (16) 
and may be rewritten here as 

q, = dE. - dE = a[dE. - dE ] . 
f i r i w 

The contribution of heat transfer to the surface by the molecules incident 
after one collision may be expressed as 

• (t ■ f) (dE r - < e ;>- (t - ?) < dE r • d V a - (40) 


Successive contributions of the heat transfer for multi-collisions 
greater than one are 


\ o\ 

3 " *) 


(l/2) j_1 





(1/2) J -1 



dE ) 
w 


(41) 
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where j = 2,3,4, . .., j Q , and represents the number of interreflections. 

The superscript numbers are used to define the energy carried away by 
the molecules caused by the j* 1 * 1 interreflection. Desired is the total 
heat transfer to a differential surface element caused by an infinite 
number of interreflections, or a finite number that fully represent the 
effects of concavity as a function of dE^ and dE w . By assuming a finite 
number of interreflections, the resulting simultaneous equations can be 
solved, producing the heat transfer equation for the incident and each 
interref lected flux of molecules. A summation of these contributions 
produce the total heat transfer to the surface element, and may be expressed 
as 



j=l 


(42) 


where q^ is defined by equation (19), and q x is obtained from equations 
(16) and (43). 


qi 



(43) 


The additional contributions of the heat transfer for integer values 
of j i 2 is similarly obtained and may be written 


= q £ (l - U/2)J- 1 (I - f) 

where j = 2, 3, 4, j . The total heat transfer to a differential 

surface element in region I is obtained by inserting equations (43) and 
(44) into equation (42) 
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(45) 




j=l 


(1 - a) j 


( 1 / 2 ) 



The heat transfer to a differential surface element in region II is 
obtained similarly by excluding the freestream contribution. 

- /> (1 . a)J (1/2) J_i . (46) 

f=i 



The average heat transfer to the surface is thus obtained by integrating 
equations (45) and (46) over the surface 


Q = a 


q I dA I A 


+ i r 


j 

n 


q n dA ir 


(47) 


Values of the accommodation coefficient determined experimentally for air 
on various surfaces are shown to be between 0.88 and 0.94 [4], It is 
therefore considered justified to assume that a = 1. For this assumption, 
the average heat transfer equation reduces to 



(48) 


where A is the surface area, and is defined by equation (19). By 
inserting these relationships into equation (48), the average heat 
transfer equation becomes 
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jt-a 


Q = — n/rT/2j 

x it 


S 2 + 


- 1 



f (e) de 


+&*** r 

a 


-S 2 s in 2 0 
e 


de 


( 49 ) 


where y is the ratio of specific heats and a function of the type of 
gas (for air, / = 7/5), and f ( 0 ) is defined by equation (34). For con- 
venience, equation (49) may be made dimensionless by dividing by 
p(2RT) 3 / 2 . The first equation for the heat transfer coefficient may 
then be written 


C 

q 


4^ 




y + 1 
2( 7 - 1 ) 


rp tt-OL 

Zsi r 

Tj J 

a 


f (e) de 


+ 


8^ 


jt-a 



a 


(50) 


VI. CONCLUSIONS 

The analysis performed here can be used to determine the aerodynamic 
forces and average heat transfer to a concave hemispherical surface in 
free molecule flow. The resulting equations, although not in closed 
form, can be integrated numerically using standard computer techniques. 

As shown in the analysis, the drag and lift coefficients are good for 
an infinite number of collisions where the average heat transfer equa- 
tion is independent of the number of interreflections. 
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The resulting equations are general in that they are a function 
of surface velocity, ambient and wall temperatures, type of gas , and 
angle of attack. The equations are limited to high molecular speed 
ratios where molecules impinge only on the concave surface. 

The drag, lift, and heat transfer coefficients have been obtained 
by solving the resulting equations numerically for some typical values 
of the molecular speed ratio and surface- to-ambient- temperature ratio 
and plotted as a function of angle of attack (figures 6 through 14). 

As shown, the average heat transfer coefficient is independent of the 
concavity. Maximum values of the heat transfer and drag occur at zero 
angle of attack, decreasing with increasing angle of attack to values 
of zero at ninety degrees angle of attack. By comparison, the drag for 
the concave hemisphere is shown to be slightly higher (about one per 
cent) than the drag for a convex hemisphere. The total lift is shown 
to be very small and would be zero except for interref lected molecules. 
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Figure 1 - Cartesian Coordinates of a Differential 
Surface Element 




Figure 2 - Illustration of Lambert's Cosine Law 
of Diffusion 



Figure 3 - Diffuse Emission Between Two Surface Elements 
of Arbitrary Orientation 
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Figure U - Angular and Regional Definitions of Hemisphere 
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Figure 5 - Diffuse Emission Inside Hemisphere 
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Figure 7. Drag Coefficient, T^/t =1.0 



Figure 8. Drag Coefficient, 
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Figure 11. Lift Coefficient 
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Figure 
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Figure 13. Average Heat Transfer Coefficient. 
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Figure 14. 
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